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METHOD AND APPARATUS FOR ANALYZING DEFECT 
DATA AND A REVIEW SYSTEM 



BACKGROUND OF THE INVENTION 

The present invention relates to a defect 
data analyzing method and apparatus and a review system 
for analyzing defect distribution state from defect 
5 data detected by an inspection apparatus in a 
semiconductor wafer manufacturing process. 

In a semiconductor device manufacturing 
procedure forming a circuit pattern on a semiconductor 
substrate (hereinafter, referred to a semiconductor 

10 wafer manufacturing procedure) , a pattern defect 

inspection or a foreign matter inspection is performed 
after each process and inspection result is analyzed, 
so as to improve the yield and stability. As one of 
such methods, a defect distribution state analysis is 

15 known. Most of defect distribution deflection is 

caused by a device failure or a process error. For 
this, it is tried to identify a defect cause by 
identifying a defect distribution pattern 
characteristic to a device failure or a process error. 

20 There are conventional techniques concerning 

this as follows. JP-A-10-214866 discloses a method for 
recognizing a region where defects are concentrated in 
the defect distribution as a cluster and selecting a 
point of reviewing (viewing the detected defect by 
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using an optical microscope or a scanning electron 
microscope) according to the area and shape of the 
cluster. JP-A-6-61314 discloses a method for dividing 
wafers into groups according to the state how the 
5 defect map is clustered and judging whether the state 
is similar to a known pattern so as to identify a 
defect cause. Moreover, JP-A-11-45919 discloses a 
method for identifying a defect cause by matching 
defect distribution image data with case database whose 

10 defect cause can be identified. U.S. Patent No. 

5,982,920 discloses a method for classifying defects 
into user-defined events associated with defect causes 
according to the defect distribution state. 

According to the method disclosed in the 

15 aforementioned JP-A-11-45919 , defect distribution image 
data is created by using the pixel value of the image 
data corresponding to the wafer region as a value 
proportional to the number of defects in the region 
corresponding to pixels. From the defect distribution 

20 image data, a blob of defect concentration, i.e., a 
cluster is detected. The cluster portion is matched 
with the case database. According to the information 
associated with the data of a high matching degree, a 
defect cause is identified. 

25 According to the method disclosed in U.S. 

Patent No. 5,982,920, according to the shape of the 
defect-concentrated portions, the concentrations of 
defects are classified into microstructure clusters, 
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curvilinear clusters, and amorphous clusters and the 
other portions are classified as global. They. are 
respectively classified into user-defined events and 
related to defect causes. 
5 In the aforementioned methods, it is possible 

to recognize a cluster when a remarkable defect 
distribution pattern appears. However, in JP-A-11- 
45919, no consideration is taken on a case when the 
defect distribution pattern to be recognized is weak, 

10 i.e., when the defect density is low or when the 

difference between the inside and outside of the patter 
is small. According to the method disclosed in U.S. 
Patent No. 5,982,920, such cases are classified as 
global. However, various shapes are involved there and 

15 should be classified into user-defined events so as to 
identify a defect cause. 



SUMMARY OF THE INVENTION 

It is therefore an object of the present 
invention to provide a defect data analysis method and 

20 device and a review system capable of coping with a 
weak defect distribution pattern and outputting 
information supporting to identify a failure cause at a 
high speed without requiring user definitions. 

That is, the present invention provides a 

25 method for analyzing defect distribution from defect 

data obtained by inspecting a substrate processed in a 
process for forming a circuit pattern on a substrate. 



From the information on the defect position obtained by 
inspecting the processed substrate, distribution of the 
defect on the processed substrate is obtained and by 
using information on defect positions on the processed 
substrate, characteristic of the obtained defect 
distribution is classified into one of the following 
distribution characteristic categories: repeated 
defects, clustered defects, arc-shaped regional 
defects, radial regional defects, line type regional 
defects, ring and blob type regional defects, and 
random defects. 

These and other objects, features and 
advantages of the invention will be apparent from the 
following more particular description of the 
embodiments of the invention, as illustrated in the 
accompanying drawings . 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing a concept of a 
defect data analysis method according to the present 
invention . 

Fig. 2 is a flow for identifying repeated 

defects . 

Fig. 3 is an example of the nearest-point 
Voronoi diagram. 

Fig. 4A and Fig. 4B explain the principle of 
detecting an arc. 

Fig. 5 is a flow for identifying arc-shaped 
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regional defects. 

Fig. 6 explains the principle of the Hough 

transform. 

Fig. 7 shows an example of combinatorial 
5 Hough transform processing. 

Fig. 8 is a flow for identifying radial 
regional defects. 

Fig. 9 is a flow for identifying line type 
regional defects. 
10 Fig. 10A to Fig. 10F explain the method for 

identifying the line type regional defects. 

Fig. 11 is a diagram showing a concept of 
identification of ring and blob type regional defects. 

Fig. 12 a flow for creating a high-density 
15 region image. 

Fig. 13 explains automatic creation of a 
ring-shaped template pattern. 

Fig. 14A to Fig. 14G are diagrams explaining 
the method for automatic creation of a blob-shaped 
20 template pattern. 

Fig. 15 shows diagrams for explaining the 
method for calculating matching score of the high- 
density region image with the template pattern image. 

Fig. 16 shows diagrams for explaining the 
25 concept of ring and blob type distribution applicable 
to a partially inspected wafer. 

Fig. 17 explains a first user pattern 
registration method . 
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Fig. 18 explains a second user pattern 
registration method . 

Fig. 19 is a first configuration example of 
the inspection device according to the present 
5 invention . 

Fig. 20 is a second configuration example of 
the inspection device according to the present 
invention. 

Fig. 21 shows an example of a defect data 
10 analysis result display screen . 

Fig. 22 shows a configuration example of the 
review system according to the present invention. 

DESCRIPTION OF THE EMBODIMENTS 

Description will now be directed to 
15 embodiments of the present invention with reference to 

the attached drawings. 

Fig. 1 is a diagram showing a concept of a 

defect data analysis method according to the present 

invention . 

20 In the defect data analysis method according 

to a first embodiment of the present invention, defect 
data 102 output from a semiconductor wafer inspection 
apparatus 101 includes at least position coordinates of 
defects. 103a and 103b are wafer maps showing the 

25 defect positions on a coordinate system having one 

point on the waver as an origin. The defect position 
coordinates output from the inspection apparatus may be 
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described by a coordinate system having one point on 
the wafer as an origin or a coordinate system having an 
origin for each chip. In the former case, X and Y 
coordinate values are directly used while in the latter 
5 case, the position coordinates in the chip should be 
converted into another coordinate system by using the 
chip arrangement information and the chip size 
information. According to the defect distribution 
state on the wafer map, defects are classified into the 

10 following defect distribution characteristic 

categories: repeated defects, clustered defects, arc- 
shaped regional defects, radial regional defects, line 
type regional defects, ring and blob type regional 
defects, and random defects. 

15 The repeated defects are defects which are 

located at almost identical positions within a chip 
over a plurality of chips. The repeated defects are 
considered to be caused by reticle defects or fault 
alarm. The black points on 104a and 104b represent 

20 repeated defects corresponding to 103a and 103b, 
respectively . 

Fig. 2 shows a flow for identifying the 
repeated defects. Firstly, from the defect position 
coordinates, a chip-superimposed map is created by 

25 representing all the defect positions by chip origin- 
referenced coordinates (S201). When the defect 
position coordinates are described by an entire wafer 
coordinate system, it is necessary to perform 



conversion into the chip origin coordinate system by 
using the chip arrangement information and the chip 
size information. Next, a nearest-point Voronoi 
diagram is created to obtain a distance between 
adjacent defects on the chip-superimposed map (S 202) . 

Fig. 3 shows an example of the nearest-point 
Voronoi diagram. The nearest-point Voronoi diagram 
shows a region of influence of each coordinate point. 
The region of influence of coordinate point P is a set 
of points whose nearest coordinate points are P. This 
is expressed by a polygon surrounding the P, which is 
called a Voronoi cell of P. a point corresponding to 
the Voronoi cell adjacent to the Voronoi cell of P is 
an adjacent point of P. Next, clustered defects, i.e., 
adjacent defects in the same chip present at a distance 
d not greater than a predetermined threshold value dl 
are excluded from the object to be processed (S203) . 
Among the remaining defects, clustered defects, i.e., 
defects present at a distance d not greater than a 
predetermined threshold value d2 are combined into a 
group and the group consisting of a predetermined 
number or more of defects is detected as a repeated 
defect (S204). The repeated defects are numbered for 
each group. The in-chip coordinates, the number of 
constituting defects, and the number of chips 
characterizing the group are calculated and recorded. 

Clustered defects are characterized in that 
the distance between the adjacent defects on the wafer 
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map is very small. The clustered defects are 
considered to be caused by defocus and micro-scratch. 
The portions indicated by the black points 105a and 
105b are clustered defects corresponding to 103a and 
5 103b. In order to identify the clustered defects, 

defects at a distance not greater than a predetermined 
threshold value are connected into a group and a group 
consisting of a predetermined number of defects is 
detected as a clustered defect. A nearest-point 

10 Voronoi diagram corresponding to a wafer map is created 
to obtain a distance between adjacent defects, thereby 
checking whether the distance between the adjacent 
defects is equal to or smaller than a threshold value. 
Alternatively, it is possible to divide the wafer map 

15 into grids and check to which grid each point belongs, 
thereby checking whether the distance between the 
defects in the same grid or adjacent grids is equal to 
or smaller than the threshold value. Clustered defects 
are numbered for each group and the following factors 

20 characterizing the group are calculated and recorded: 
coordinates of center of gravity, maximum x/y 
coordinates, minimum x/y coordinates, area, defect 
density, number of defects, and so on. 

Arc-shaped regional defects are defects 

25 distributed in the arc shape due to a scratch in the 

CMP (chemical mechanical polishing) process. Since the 
wafer maps 103a and 103b have no arc-shaped regional 
defects, an example of another wafer 109 is shown. 



Radial regional defects are stripe-shaped 
defects extending from the center of the wafer toward 
outside. The radial regional defects are caused by 
foreign matters attached to a blow-out opening of a 
device for blowing gas from the center. Since the 
wafer maps 103a and 103b have no radial regional 
defects, another example of wafer 110 is given. 

Line type regional defects are defects having 
the highest defect density distributed in a straight 
line. The line type regional defects are considered to 
be caused by a convey scratch or flow. Rectangular 
portions of 106a and 106b are line type regional 
defects corresponding to 103a and 103b. The line type 
regional defects are not arranged strictly in a line 
but have a certain width. 

Ring and blob type regional defects are 
defects having a portion of the higher defect density 
distributed in a ring or blob shape. The ring and blob 
type regional defects are considered to be caused by a 
periphery peeling off or device stain. The portions 
surrounded by solid lines in 107a and 107b are ring and 
blob type regional defects corresponding to 103a and 
103b. 

Defects which have not been classified into 
any of the aforementioned defects are extracted as 
random defects. 108a and 108b show random defects 
corresponding to 103a and 103b. 

Next, detailed explanation will be given on 



1 •■ 

-li- 
the method for identifying the arc-shaped regional 
defects according to the present invention with 
reference to Figs. 4A and 4B and Fig. 5. 

Fig. 4A and Fig. 4B explain the principle of 
5 the arc-shaped defect detection. As shown in Fig. 4A, 
when arbitrary two points on the arc shape are selected 
and a perpendicular bisector is drawn, the 
perpendicular bisectors intersect one another at one 
point, which is the center of the arc shape. From 

10 this, by voting points where the perpendicular 

bisectors of all the combinations of two points from 
the point group on the xy space pass to (x, y) , it is 
possible to detect the maximum voting value (xO, yO) as 
the candidate of center of the arc. At this point of 

15 time, information on the radius r is lost and all the 

points are polar-coordinate-converted into (r, 9) using 
the center candidate (xO, yO) as the origin as shown in 
Fig. 4B. Since the radius r is constant, the point 
string on the arc in the xy space is distributed on a 

20 horizontal line in a polar coordinate space. 

Accordingly, by detecting a horizontal segment 
corresponding to an arc and inversely converting it 
into the xy space, it is possible to detect an arc. 

Fig. 5 is a flow for identifying the arc- 

25 shaped regional defects according to the aforementioned 
principle of the present invention. 

Firstly, two points are selected from the 
defect coordinates and voting is performed to a point 
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where the perpendicular bisector of the points pass 
(S501) . The voting space is twice or thrice larger 
than the wafer map. The selection of two points may 
include all the combinations but in order to reduce the 
5 time, they can be thinned at random. An appropriate 
pre-processing may be performed to exclude points 
judged to be not constituting the arc. For example, 
the aforementioned repeated defects are excluded. 
Moreover, among the clustered defects, defects of a 

10 blob having a very large width as compared to the arc 
shaped distribution is excluded. When performing 
voting, points too near from each other cause a large 
error and points too far from each other may have no 
relationship. Accordingly, weighting is performed 

15 according to the distance between the two points so as 
to reduce the contribution of the points too near or 
too far from each other. 

Next, the maximum voting value (x, y) is 
detected as a center candidate (S502) . Peak sharpness 

20 is checked by comparing with the surrounding values and 
if the peak is not sharp, the processing is terminated 
(S503) . If the peak if sharp, by using the maximum 
voting value (x, y) as the origin, all the defect 
coordinates are converted into (r, 0) so as to create a 

25 polar coordinates converting image showing the defect 
density on the r9 space as a gray level (S504). 
Projection onto the r axis is calculated and the r of 
the maximum projection value is found (S505) . 



- 13 - 

Next, a threshold value is calculated by an 
appropriate method such as a judgment and analysis 
method and a percentile method and the polar 
coordinates converting image is digitized (S506) . At 
5 the maximum voting value r, a horizontal segment is 
detected (S507). Discontinuity not greater than a 
predetermined threshold value is assumed to be 
continuous and a continuous section not smaller than a 
threshold value is assumed to be a segment. If no 

10 segment is detected, the processing is terminated 

(S508). The segment detected is inversely converted 
into the xy space to obtain an arc. Defects present at 
a distance within a threshold value from the arc are 
detected as arc-shaped regional defects (S509) . The 

15 arc-shaped regional defects detected are excluded from 
the object to be processed (S510) and control is 
returned to S501 to repeat the processing. 

In this embodiment, explanation has been 
given on the method for detecting the center by voting 

20 perpendicular bisectors to the xy space. However, 
other methods ca be used to detect the center. For 
example, there are a method for voting traces of 
circles passing through a single point to a parameter 
space of circles consisting of (x, y, r) , and a method 

25 for selecting arbitrary three points, calculating a 

circle passing through these points, and voting it to 
the (x, y, r) parameter space. 

Next, detailed explanation will be given on a 
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method for identifying radial regional defects with 
reference to Fig. 6 and Fig. 7. 

For identification of the radial regional 
defects, it is possible to use a combinatorial Hough 
5 transform which is an improved Hough transform. 
Firstly, explanation will be given on the Hough 
transform. Fig. 6 explains the principle of the Hough 
transform. A straight line on the xy space is 
expressed by two parameters: a distance c from the 

10 origin and the angle 0 of the perpendicular line. When 
a straight line passing through point PI on the xy 
space is plotted onto the p9 space, the trace is a curve 
passing through (91, pi) and (92, p2 ) as shown in the 
figure. For the other points on the xy space, it is 

15 possible to draw curves in the same way. When points 
are distributed on a straight line as shown in the 
figure, their traces intersect one another at (92, p2) 
representing a straight line passing through these 
points. Accordingly, p9 space traces of points on the 

20 xy space are voted onto the (9, p) and a straight line 
corresponding to the maximum voting number (9, p) is 
detected . 

The combinatorial Hough transform is also 
voting to the p9 space but voting is performed only to 
25 one point on the p9 space corresponding to a stright 
line passing through two points on the xy space. If 
the coordinates of the two points PI and P2 on the xy 
space are (xl, yl) and (x2, y2 ) , then the corresponding 
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p and 0 can be calculated by the following expressions. 



0 = 



tan -1 ((xl - x2)/(y2 - yl))+ it (if y2 - yl < 0 D xl - x2 > 0 ) 
tan- 1 ((xl-x2)/(y2-yl))-7c (if y2 - yl < Of) xl - x2 < 0 ) but |xl|>|x2| 
tan" {{x\-x2)l(y2-y\)) (eles) 



(Expression 1 ) 



p = cos<9 • (xl + x2)/2 + sin6> • (y\ + y2)/2 (Expression 2) 

Here, PI is the one having a greater absolute 
value of the x coordinate. In the normal Hough 
transform, the is in the range from -n/2 to tt/2 but here 
5 it is in the range from -71 to n. That is, hot only the 
inclination of the straight line but its direction is 
also considered. Fig. 7 shows correspondence between 
combinations of two points on the xy space and the 
points on the p0 space. Like symbols correspond to each 
10 other. As shown in the figure, when radial regional 

defects are present, a plenty of peaks are detected at 
different 0 in the vicinity ofp = 0. This is utilized 
for identification of the radial regional defects. 

Fig. 8 is a flow for identifying radio 
15 regional defects according to the present invention. 

Firstly, arbitrary two points are selected on 
the defect coordinates and voting is performed to the 
point on the p6 space corresponding to a straight line 
passing through the two points (S801) . The method for 
20 selecting the two points may be all the combinations 



but in order to reduce the time, they can be thinned at 
random. Moreover, points in the vicinity of the center 
of the wafer are inevitably voted to the vicinity of p 
= 0 and become a noise for the voting from the radial 
regional defects. Accordingly, it is possible to 
exclude defects whose distances from the center are 
equal to or below a predetermined value. When 
performing voting, in the same way as the 
identification of the arc-shaped regional defects, the 
distances between two points are weighted so as to 
reduce the contribution of the points too near or too 
far from each other. 

Next, a predetermined number of peaks are 
detected in the descending order of the values (S802) . 
When arc-shaped regional defects are present, voting is 
concentrated around p = 0. Accordingly, in order to 
check the concentration degree, the ratio a of peaks in 
a predetermined range around p = 0 against all the 
number of detected peaks is calculated (S803) . The a 
may also be a ratio of voting performed to a 
predetermined range around c against all the voting. 
When the a is equal to or below a predetermined value, 
it is judged that no radial regional defects are 
present and the processing is terminated (S804) . 

Next, an appropriate threshold value is 
decided for the peaks present within a predetermined 
range around p = 0 and the vicinity of the peaks, for 
example, an area of only 5 x 5 is digitized and labeled 



(S805) . The threshold value may be decided to be a 
constant ratio against the peak value or decided from 
the higher value of the aforementioned vicinity area to 
make a predetermined number of pixels "1". If the 
5 number of labels is equal to or above a predetermined 
number, it is decided that radial regional defects are 
present (S806) . The ranges of p and 0 are decided from 
each label size and corresponding defects on the xy 
space are detected as radial regional defects (S807). 
. 10 Next, detailed explanation will be given on 

identification of line type regional defects with 
reference to Fig. 9 and Fig. 10A to Fig. 10F. 

Identification of line type regional defects 
according to the present invention uses the Hough 

15 transform. Fig. 9 shows a flow for identifying line 
type regional defects. Fig. 10A to Fig. 10F show 
examples of processing. Firstly, defect coordinates 
are subjected to the Hough transform (S901) and (0, p) 
of the maximum voting count is detected as a candidate 

20 straight line (S902) . Fig. 10B is a Hough transform 

image corresponding to the wafer map shown in Fig. 10A. 
As has been described above, the line type defects have 
a certain width and the resolution of (0, p) is made 
rough. However, if normal Hough transform is 

25 performed, for example, there arises a problem that 

when clustered defects are present at two positions, a 
straight line connecting the two clustered defects is 
detected. The clustered defects may be distributed in 
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a linear shape and accordingly, exclusion of the 
clustered defects adversely affects the result. 

In order to solve this problem, contribution 
of the clustered defects is reduced by performing 
5 weighting proportional to the distance between defects 
or square of the distance between the defects when 
voting. For this, the nearest-point Voronoi diagram 
corresponding to the wafer map is created in advance 
and the distance between the defects is calculated. If 

10 a Voronoi diagram has been created when identifying the 
clustered defects, it can be used. Moreover, since one 
defect is present in a Voronoi cell, a reciprocal 
number of the area of the Voronoi cell may be 
considered as a local defect density in the 

15 corresponding coordinates. For the aforementioned 

weighting, it is possible to use the reciprocal number 
of the defect density, i.e., the Voronoi cell area. 

Next, as shown in Fig. 10B, a width W of the 
pdirection of a portion not smaller than a % of the 

20 maximum voting count around the candidate straight line 
is calculated (S903) . The a is a predetermined 
threshold value and W is considered to be a width of 
the candidate straight line. Next, check is made to 
determine whether the W is not greater than the 

25 predetermined threshold value. When the width is 

greater than the threshold value, i.e., when the width 
is too great, it is judged not line type regional 
defects are present and the processing is terminated 



(S904). When the width is smaller than the threshold 
value, the wafer map is rotated by 9 and a gray level 
image representing the density around the candidate 
straight line is created (S905) . Here, the image size 
5 is made identical to W. 

The grating shown in Fig. IOC corresponds to 
the pixels of a gray level image. It should be noted 
that in the figure, instead of rotating the wafer map 
by -0, the grating is rotated by 0. The gray level 

10 image is created by corresponding a greater pixel value 
to the grid having a high defect density and a small 
pixel value to the grid having a low defect density. 
For example, the pixel value is proportional to the 
number of defects in the grid. It is possible to 

15 perform weighting such that performed upon the Hough 
transform processing in order to reduce the 
contribution of the clustered defects. 

Next, by using an appropriate method such as 
judgement and analysis method, the gray level image is 

20 digitized as shown in Fig. 10D and from the digitized 

image, it is decided whether line type regional defects 
are present (S906) . The portions shown in black in 
Fig. 10D are candidate of the line type regional 
defects. Since the center row corresponds to the 

25 position of the candidate straight line, the 

presence/absence of the line type regional defects is 
decided by the length of the candidate portion of the 
central row, whether the candidate portion is 
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continuous, by presence/absence of distribution in the 
right and left rows converted into numeric values which 
are compared to a predetermined threshold value. 

In the example shown in Fig. 10D, the 
5 rectangular portion indicated by a solid line in Fig. 
10E can be judged to be a line type distribution. When 
no line type distribution is present, the processing is 
terminated (S907) . When the line type distribution is 
present, as shown in Fig. 10F, Fig. 10E is rotated by 0 

10 so as to be superimposed on the corresponding position 
of the wafer map and the defects inside the solid 
rectangular shape are detected as line type regional 
defects (S908) . 

The line type regional defects detected are 

15 excluded from the object to be processed (S909) and 

control is returned to S901 to repeat the processing. 
For each line type regional defects which are detected 
by one processing is added by a group number and its 
position, width, angle, length, defect density, and 

20 other characteristic amounts are calculated and 
recorded . 

Next, detailed explanation will be given on 
the method for identifying the ring and blob type 
regional defects according to the present invention. 
25 Fig. 11 shows a concept of identification of the ring 

and blob type regional defects according to the present 
invention. In the first step, a digital image is 
created by expressing higher defect density portions by 
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1 and lower defect density portions by 0 according to 
the defect position coordinates. Hereinafter, this 
image will be referred to as a high density area image 
704. In the next step, a plurality of geometric 
5 template patterns 705 are automatically created in 

accordance with the image size of the high density area 
image 704. Next, matching scores are calculated 
between the high density area image 704 and the 
automatically created plurality of geometric template 
10 pattern images 705. In the last step, a pattern image 
706 having the highest matching score is selected. 

According to another embodiment of the 
present invention, in addition to the automatically 
created geometric template patterns 705, a user 
15 registers a template pattern 707 in advance and a 

pattern image 706 having the highest matching score is 
selected. The pattern image 706 is superimposed on the 
wafer map 103 and defects contained in the pattern 
portion are detected as ring and blob type defects. 
20 Next, each step will be detailed. 

Firstly, explanation will be given on the 
method for creating the high density area image 704 
according to the defect position coordinates. 

Fig. 12 shows a flow for creating the high 
25 density area image 704. 

Firstly, a nearest-point Voronoi diagram is 
created for coordinates of all the defects and an area 
of Voronoi cell of each point is calculated (S1201) . 
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If this calculation has been performed in the 
aforementioned line type regional defect 

identification, the results can be used here. Next, a 
histogram of the area of the Voronoi cell is calculated 
5 (S1202). It should be noted that the defects which 
have been classified into the repeated defects, 
clustered defects, or the line type regional defects 
are not included in the histogram. A modal value is 
obtained from the histogram and the modal value is 

10 multiplied by a predetermined number (N) . Its square 
root is calculated to be a pixel size (S1203). The 
image size is determined so that the entire wafer is 
contained just in one image and an image where the 
number of defects per pixel is represented by the gray 

15 level is created (S1204). Such a image can be created 
by initializing all the pixels to 0 and incrementing 
the pixel value of the position corresponding to the 
coordinates of each defect. Here, if the N is too 
large, the image size is small and it becomes difficult 

20 to recognize the shape. The N is preferably about 5. 
On the contrary, if the N is too small, the difference 
between the gray levels between the high density 
portion and the low density portion becomes small and 
it becomes difficult to recognize the difference. 

25 Next, a predetermined threshold value (T) is 

used to obtain a digitized image in which high density 
portions are expressed by 1 and low density portions 
are expressed by 0 (S1205) . Here, the value T is 



preferably similar to the N. In the image obtained, a 
pattern of relatively high defect density can be seen. 
Lastly, the pattern is subjected to 
expansion/contraction processing to obtain a high 
5 density area image 1104 (S1206) . By creating an image 
by such a method, even if the distribution density is 
low, it is possible to set an appropriate image size 
and extract a relatively high density area. Moreover, 
another method can be used by creating a gray image in 

10 S1204 and performing digitization in S1205. Each 
defect is weighted proportionally to the shortest 
distance to the adjacent defect or its square or 
Voronoi cell area and it is added to the pixel value at 
the corresponding position so as to create a gray image 

15 and perform digitization by using the judgment and 
analysis method 

Next, explanation will be given on the 
automatic creation of the geometric template pattern 
image 7 05. 

20 The geometric template pattern image 705 is 

automatically created separately for the ring pattern 
and blob pattern. The image size is identical to the 
high density area image 704 and the pattern is created 
as follows by using the maximum circle that can be 

25 drawn on the image as a reference. 

The ring pattern is created by combination of 
small areas obtained by equally dividing the reference 
circle into concentric circles in the radial direction 
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and equally dividing by sectors in the center angle 
direction. Fig. 13 shows an example of pattern divided 
into five portions in the radial direction and eight 
portions in the center angle direction. Explanation 
5 will be given on the method of combination of the small 
areas using an example of a pattern divided into five 
parts in radial direction and eight parts in the center 
angle direction. Firstly, symbols a, b, c, d, and e 
are assigned to the concentric annuluses (the innermost 

10 one is a circle) obtained by dividing the circle in the 
radial direction in the order counted from the outside 
and numbers 1 to 8 are assigned to the sectors obtained 
by dividing the circle into eight equal sectors in the 
center angle direction. There are fifteen combinations 

15 of complete annulus or circle size: a, b, c, d, e, ab, 
be, cd, de, abc, bed, cde, abed, bede, and abede . As 
for the part variations of the annuluses, there are 
eight combinations for each of the sizes: one eighth, 
two eighths, three eighths, four eighths, five eighths, 

20 six eighths, seven eighths. For example, for the size 
of two eighths, there are combinations of 12, 23, 34, 
45, 56, 67, 78, 81. By adding a complete annulus to 
these combinations, 57 combinations can be considered. 
For each of them, the aforementioned 15 size 

25 combinations can be considered. Accordingly, 855 
combinations of patterns are generated in total. 

The division in the radial direction and the 
center angle direction may be different from this. 
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Moreover, the division may be performed not equally. 
However, the equal division is preferable to simplify 
the processing. 

The blob type pattern is generated by a 
5 combination of small areas obtained by dividing the 

reference circle into ellipses having major axes in the 
vertical direction equally in the horizontal direction 
and dividing the reference circle into ellipses having 
major axes in the horizontal direction equally in the 

10 vertical direction. 

Fig. 14A shows an example of dividing the 
reference circle into eight ellipses in the horizontal 
and the vertical directions. There are eight size 
combinations from 1 to 8 and (9-N) combinations of 

15 positions when the size is N in the horizontal 

direction. Accordingly, there are 36 combinations in 
total. There are also 36 combinations in the vertical 
direction. Consequently, 1296 patterns can be 
generated by combinations of the horizontal direction 

20 and the vertical direction. Fig. 14B to Fig. 14 D show 
examples of patterns generated by this method. The 
black portions correspond to the patterns. A 
combination of separate areas as shown in Fig. 14E or 
the irregular pattern as shown in Fig. 14F are not 

25 generated by this method. The blob type pattern may 

also be obtained by rotating the reference circle by 45 
degrees around the its center. Fig. 14G shows the 
areas obtained in this case. It is possible to use one 
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of the reference circles or both of them. Moreover, 
the number of divisions may be other than eight. 

By automatically generating the geometrical 
template pattern image 705 by the aforementioned 
5 methods, it is possible to reduce the number of pattern 
to a practical level, thereby realizing a high speed 
processing . 

Next, explanation will be given on the method 
for calculating the matching score between the high 

10 density area image 704 and the template pattern image 

705 with reference to Fig. 15. Each pixel value of the 
high density area image 704 is compared to the pixel 
value of the identical address of the template pattern 
image 705. If matched, +1 is obtained and if not 

15 matched, -1 is obtained. The comparison results of all 
the pixels are totaled and the obtained value is used 
as the matching score. In Fig. 15, the high density 
area image 704 and the template pattern image parts 
705a and 705b are shown by expressing the pixel value 1 

20 by black and pixel value 0 by white. 708a and 708b 

show unmatched portions between the high density area 
image 704 and the template pattern images 705a and 705b 
by black. That is, in these images, the pixel value of 
the white portion is 1 and the pixel value of the black 

25 portion is -1. The total of all the pixel values is 
the matching score. When 708a is compared to 708b, 
708a has a smaller black area and accordingly, 705a has 
a higher matching score with the high density area 



image 704 than 705b. 

According to the method explained above, the 
pattern 706 is detected without fail. Actually, 
however, when the density is very low or the density 
difference is very small, it is better not to obtain 
any pattern. Accordingly, the detected pattern 706 is 
verified to judge whether a pattern is present. 
Hereinafter, explanation will be given on the judgement 
method for presence/absence of a pattern and the 
pattern detection sensitivity adjustment method. In 
order to judge whether a pattern is present, a gray 
level image is created by the method for performing 
weighting when creating the high density area image 704 
and the gray level image obtained is used. The 
detected pattern image 706 is superimposed on the gray 
level image and a ratio C of the j udgment /analysis 
value inside and outside the pattern is calculated from 
the pixel values of the gray level image by using the 
following expression . 



wherein, 

(co p , co b ; ) number of pixels inside and outside the 
pattern 

(p p , ]i b ; ) average value of pixel values of gray level 
image inside and outside the pattern 




(Expression 3) 
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(a p/ a b ; ) standard deviation of pixel values of gray 
level image inside and outside the pattern 
The C is compared to a predetermined threshold value 
and if the C is equal to or below the threshold value, 
5 it is judged that there is no pattern. Moreover, when 
the modal value in the histogram obtained in S1202 is 
lower than a predetermined threshold value, no more 
processing is performed and it is judged that there is 
no pattern. When the former threshold value is set to 

10 a small value, it is possible to detect a patter of a 
smaller density difference. Moreover, when the latter 
threshold value is set to a small value, it is possible 
to detect a pattern of a smaller defect density. If 
these threshold values can be specified by a user, the 

15 user can adjust the sensitivity. 

By the way, the semiconductor wafer 
inspection does not always inspect the entire surface 
because of the throughput restriction. It is often the 
case that a partial inspection is performed by 

20 specifying inspection/non-inspection for each chip row. 
In this case, the pattern detection cannot be performed 
easily by using the aforementioned ring and blob type 
regional defect identification. In order to solve this 
problem, another method for identifying the ring and 

25 blob type defects will be explained with reference to 

Fig. 16. Fig. 16 shows a concept of identification of 
the ring and blob type defects applicable to defect 
data of a partial inspection wafer. When a partial 
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inspection is performed, a wafer map 103 looks as shown 
in the figure, for example. In the same way as has 
been described above, the high density area image 704 
is created and a geometrical template pattern 705 is 
5 automatically created. On the other hand, according to 
the inspection condition information and the chip 
matrix data, an inspection area image 709 is created. 
The inspection area image 709 is superimposed on the 
geometrical template pattern 705. The non-inspection 

10 area is masked and the matching score is calculated. 
That is, the template pattern image 705 and the 
inspection area image 709 of the high density area 
image 704 are simultaneously scanned. If the area is 
the non-inspection area, then 0 is obtained. If the 

15 area is the inspection area and if the pixel value of 
the high density area image 704 matches with the pixel 
value of the template pattern image 705, then +1 is 
obtained. Otherwise, -1 is obtained. These values are 
totaled to obtain the matching score. Thus, the 

20 matching score with each template pattern 705 is 
calculated and a pattern image 706 of the highest 
matching score is selected. According to this method, 
even in the partial inspection, it is possible to 
accurately calculate the matching score with the 

25 template pattern 705 and accordingly, it is possible to 
identify the ring and blob type region defects. 

Next, explanation will be given on the user 
pattern registration method. 
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Fig. 17 explains the user pattern 
registration method. Firstly, a user specifies a wafer 
map 901 having the pattern to be registered and display 
it on the screen. This specification may be omitted. 
5 In this case, nothing is displayed. Next, one of the 
division patterns 902a to 902c is selected for 
automatic creation of the geometric template pattern. 
902a is identical to the ring type pattern shown in 
Fig. 13. 902b and 902c are identical to the blob type 

10 pattern shown in Fig. 14A and Fig. 14G, respectively. 
The division pattern 902 selected is superimposed on 
the wafer map 901. Moreover, each of the small areas 
of the selected division pattern 902 is switched 
between selection and non-selection states according to 

15 the user specification and the areas selected are 
displayed with a different color. 

The user selects a position having a high 
defect density or arbitrary small areas on the wafer 
map and specifies registration after the selection. 

20 The combination of the small area selected is compared 
to the pattern automatically created. If the 
combination is not included in the pattern 
automatically created, the pattern image 903 is stored 
as the user registration pattern 707. According to 

25 this method, it is enough to encode the combination of 
selected small areas and record it as the user 
registration pattern instead of using a pattern image. 

Fig. 18 explains a user pattern registration 



method which is different from the aforementioned. A 
user specifies the wafer map 901 having the pattern to 
be registered and display it on the screen. A high 
density area 904 is created according to the 
5 aforementioned high density area image creation method 
and displayed on the screen. Matching scores are 
calculated between the high density area 904 and the 
template pattern automatically created and the user 
registration pattern already registered. When the 
10 maximum value of the matching scores is not greater 

than a predetermined threshold value, it is registered 
as the user registration pattern 707. In this case, 
the user registration pattern 707 is recorded as image 
data . 

15 As the user pattern registration method, it 

is possible to provide the aforementioned two methods 
or one of them. Moreover, it is possible to use a 
method different from the aforementioned methods. 

In the defect data analysis method according 

20 to a second embodiment of the present invention, when 
identifying ring and blob type regional defects, by 
adding information associated with the geometric 
template patter 705 automatically created and the user 
registration pattern 707, the pattern image 706 is 

25 selected in the same way as the first embodiment. 

The information added includes pattern 
position, size, shape, importance, and special matter. 
The pattern position, size, and shape of the geometric 
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template patter 705 are roughly classified according to 
a certain rule and default information is added. For 
the blob type pattern, the pattern position is 
classified, for example, upper right, top, upper left, 
5 right, center, left, lower right, bottom, and lower 

left. The size is classified, for example, into large, 
intermediate, and small. The shape is a blob. For 
example, Fig. 14B has a small size blob pattern at the 
upper left, Fig. 14C has an intermediate size blob 

10 pattern at the right, and Fig. 14D has a large size 
blob pattern at the center. As for the ring type 
pattern, it is a radius-size quarter ring pattern at 
the upper right. It should be noted that even if the 
pattern derives from a blob type pattern, the pattern 

15 is handled as a ring pattern if the pattern consists of 
only one outer row, and even if the pattern derives 
from a ring pattern, the pattern is handled as a blob 
pattern if the innermost circle is contained. 

Furthermore, the geometric template pattern 

20 705 and the user registration pattern are manually 

grouped and means is provided for collectively adding 
information, so that the same information can be added 
to the patterns which the user consider identical. 
Thus, it is possible to adjust resolution of the 

25 position and size information. The information of 

importance is input by the user. The special matter 
includes information associated directly with the 
defect causes. 
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Next, explanation will be given on an 
inspection apparatus using the defect data analysis 
method according to a first embodiment of the present 
invention. In the explanation below, repeated defects, 
5 clustered defects, line type regional defects, ring and 
blob type regional defects, and random defects are 
called regional characteristic categories. As the 
semiconductor wafer inspection apparatus, there are 
known a foreign matter inspection apparatus, an optical 

10 pattern defect inspection apparatus, SEM type pattern 
defect inspection apparatus, and the like. The 
inspection apparatus of the present invention inspects 
a semiconductor wafer by one of the known methods used 
in the aforementioned inspection apparatuses, uses the 

15 obtained defect data information to classify the 
defects into the distribution characteristic 
categories, and outputs the category information 
together with the defect data information. 

Fig. 19 shows configuration of the inspection 

20 apparatus according to the first embodiment of the 
present invention . 

Pattern information addition means 301 adds 
information associated with a pattern such as pattern 
position, size, shape, importance, and special matter 

25 to a plurality of geometric template patterns 

automatically created and the user registration 
template pattern. The pattern information addition 
means 301 consists of a section for adding default 
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information to the geometric template pattern and a 
section for manually adding information to the user 
registration template pattern. The pattern information 
addition is performed offline. The addition is 
5 performed at least once and the result is stored as a 
pattern information file 311 in storage means 302 such 
as a hard disc. 

The defect data analysis may be performed in 
both of the inline and offline modes. 

10 The inspection apparatus 303 inspects the 

wafer according to a known method and stores defect 
data 102 including at least defect position coordinates 
in the storage means 302. 

In the case of inline mode, the inspection 

15 means 303 transmits notification of the inspection end 
to the defect data analysis means 304 upon completion 
of an inspection of each wafer. The defect data 
analysis means 304 reads the defect data 102 of one 
wafer from the storage means 302. Alternatively, the 

20 defect data 102 may be passed without using the storage 
means 302 . 

In the case of offline mode, the defect data 
analysis means 304 reads defect data 102 of the wafer 
specified by an operator. 
25 The defect data analysis means 304 performs 

classification into the distribution characteristic 
categories according to the defect position coordinates 
and adds a distribution characteristic category number 
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and a group number within the category to defect data 
of each defect. When ring and blob type defects are 
detected, the pattern information file 311 is read out 
from the storage means 302. A pattern image number or 
5 pattern information is added to the defect data of the 
ring and blob type regional defects and it is stored as 
defect data 102 including the inspection and analysis 
results in the storage means 302. In the case of 
inline mode, the end of the analysis is notified to the 

10 inspection means 303. 

The result of classification into the 
distribution characteristic categories is shown on the 
display means 305. Different categories may be 
displayed by different colors on the wafer map or as 

15 shown in Fig. 1, a wafer map may be created for each of 
the categories. Simultaneously with this, the chip may 
be superimposed on the map. Moreover, the pattern 
image 706 selected in identification of the ring and 
blob type regional defects is displayed together with 

20 the information associated with the wafer to be 

inspected and the pattern information added to the 
pattern on the display means 305. Simultaneously with 
this, it is related to the wafer to be inspected and 
stored as an image file. Moreover, when a line type 

25 regional defect is detected, the rectangular shape 

showing the position is superimposed on the wafer map 
and information on the position, width, angle, length, 
defect density and the like is displayed. These 



information are related to the wafer to be inspected 
and stored. 

The user pattern registration means 306 
creates a user registration pattern by an instruction 
by an operator and stores the pattern code or pattern 
image in the storage means 302 according to the 
aforementioned method difference. It is also possible 
to configure the apparatus without any user pattern 
registration means 306. 

Next, explanation will be given on an 
inspection apparatus having the defect data analysis 
method according to a second embodiment of the present 
invention. The second embodiment inspects a 
semiconductor wafer by a known method, uses the 
obtained defect data information to classify defects 
into the distribution characteristic categories 
according to the aforementioned method, performs 
sampling according to a rule specified for each 
category, and outputs the defect data, the 
classification result, and sampling result. 

Fig. 20 shows configuration of the inspection 
apparatus according to the second embodiment of the 
present invention . 

The pattern information addition means 301 
adds a sampling condition number in addition to the 
aforementioned information. By manually grouping the 
template patterns using the same sampling condition, 
the same sampling number is added. When no sampling 
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condition is present corresponding to the number, a 
sampling condition is created by sampling condition 
creation means 307. 

The sampling condition creation means 307 
creates a sampling condition to decide a defect for 
performing review. The sampling condition related to 
the ring and blob type regional defect pattern and 
sampling conditions of other distribution 
characteristic categories are manually created and 
stored as a sampling condition file 312 in the storage 
means 302. The sampling condition may be, for example, 
the sampling count or the sampling ratio and sampling 
method. The sampling count or the sampling ratio is 
selected before entering a numeric. The sampling count 
is decided by the sampling ratio for the number of 
defects in the pattern. As for the sampling method, a 
random method, a defect number sequence method and 
other considerable methods are listed up, so that the 
user can select one from them. When the defect number 
sequence is specified, sampling is performed at the 
interval in accordance with the sampling ratio. 

The inspection means 303 and the defect data 
analysis means 304 operate in the same way as in the 
first embodiment. 

Upon completion of defect analysis, the 
sampling means 308 reads in the sampling condition file 
312 and defect data 102 including the distribution 
characteristic information of one wafer from the 
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storage means 302. The defect data 102 may be passed 
without using the storage means 302. Next, operation 
is performed according to the distribution 
characteristic category to which the defect data 
belongs. IF the defects are ring and blob type 
regional defects, sampling is performed by using the 
sampling condition related to the pattern information 
added. Otherwise, sampling is performed according to 
the sampling condition of the distribution 
characteristic category. As a result of sampling, a 
flag indicating review/non-review to be performed is 
added to the defect data and stored as a defect data 
102 including the results of inspection, analysis, and 
sampling in the storage means 302. 

The pattern image 706 selected by the defect 
data analysis means 304 is displayed by the display 
means 305 together with the information related to the 
wafer to be inspected and the patter information added 
to the pattern. The defect sampled is displayed on the 
wafer map with a color or symbol different from the 
other defects. Fig. 21 shows an example of result 
display . 

The user pattern registration means 306 
operates in the same way as in the first embodiment. 
However, the apparatus may also have configuration 
without this. 

Next, explanation will be given on a review 
system using the defect data analysis method of the 
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present invention . 

Fig. 22 explains the review system according 
to the present invention . 

The review system includes a foreign matter 
5 inspection apparatus 401, an optical pattern defect 
inspection apparatus 402, an SEM type pattern defect 
inspection apparatus 403, a defect data analysis 
apparatus 404, and a review apparatus 405 connected to 
a network. The inspection apparatuses 401 to 403 and 
10 the review apparatus 405 are arranged in a clean room, 
the defect data analysis apparatus 404 may be in any 
place . 

The results inspected by the inspection 
apparatuses 401 to 403 are file-output as defect data 

15 102 in the same format and transferred to the defect 

data analysis apparatus 404. The defect data analysis 
apparatus 404 has configuration of the inspection 
apparatus shown in Fig. 20 excluding the inspection 
means 303. The defect data analysis apparatus 404 

20 reads in the defect data 102 of the wafer to be 

analyzed, classifies the defects into a distribution 
characteristic category and performs sampling according 
to the classification result. The distribution 
characteristic information and the sampling flag are 

25 added to the defect data 102, which is file-output and 
transferred as defect data 102' to the review apparatus 
405. 

The review apparatus 405 reads the defect 



data 102 1 having the sampling flag corresponding to the 
wafer specified to be reviewed, reviews the defects 
according to the information, and manually or 
automatically classifies the defects. The category 
5 number as a classification result is added to the 

defect data 102 f , a file is output, and transferred to 
the defect data analysis apparatus 404. Simultaneously 
with this, associated review image is also transferred. 

According to the defect distribution pattern 

10 and the defect review result of the wafer to be 
analyzed, the defect data analysis apparatus 404 
creates and stores a report. Information described in 
the report includes the information on the wafer as an 
object, the distribution characteristic category 

15 classification result, and the review image. Since the 
distribution characteristic and the defect review image 
are simultaneously displayed, it becomes easier to 
identify a defect cause. According to a user request, 
more detailed information may also be described 

20 including the sampling position information, defect 
type ratio calculated from the review result and the 
sampling ratio, and other information obtained by known 
analysis means. The report file format is not limited 
to a specific one. The report file may be formatted in 

25 the HTML format which is uploaded in a data server 406 
connected by the Internet or Intranet so that it can be 
read from an arbitrary terminal 407. 

The review system which has been explained 
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above includes three types of inspection apparatuses 
and one review apparatus. However, the system may 
include only one inspection apparatus and one review 
apparatus. Moreover, the system may include a 
5 plurality of inspection apparatuses and review 
apparatuses. When the inspection apparatus has 
configuration identical to the inspection apparatus 
according to the second embodiment of the present 
invention, the function of the defect data analysis 

10 apparatus can be realized in the inspection apparatus 
and it is possible to exclude an independent defect 
data analysis apparatus. 

The aforementioned defect data analysis 
method, the inspection apparatus or the review system 

15 using the method assume that defect data of one wafer 
is to be processed. However, when not sampling but 
identification of defect distribution is to be 
performed, a plurality of wafers may also be analyzed. 
In this case, defect position coordinates of the 

20 plurality of wafers are superimposed and a high density 
area image is created according to the superimposed 
defect position coordinates. 

According to the present invention, it is 
possible to easily identify a defect cause since 

25 according to the defect position coordinates inspected 
by the inspection apparatus, the defects are classified 
into distribution characteristic categories attributed 
to different generation causes: repeated defects, 



clustered defects, arc-shaped regional defects, radial 
regional defects, line type regional defects, ring and 
blob type regional defects, and random defects. 

Moreover, since a high density area image 
representing the high defect density portion is created 
with an optimal size, it is possible to extract a 
pattern area even if the pattern is a very weak defect 
distribution pattern . 

Moreover, according to the present invention, 
the high density area image is classified into one of 
the plurality of geometric patterns and accordingly, it 
is possible to realize a defect distribution state 
analysis at a practical level of the calculation time 
and the storage capacity. 

Furthermore, since it is possible to define 
each of the plurality of geometric patterns, it is 
possible to easily identify the defect cause including 
the device and process failure. 

The invention may be embodied in other 
specific forms without departing from the spirit or 
essential characteristics thereof. The present 
embodiment is therefore to be considered in all 
respects as illustrative and not restrictive, the cope 
of the invention being indicated by the appended claims 
rather than by the foregoing description and all 
changes which come within the meaning and range of 
equivalency of the claims are therefore intended to be 
embraced therein. 



